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VERY  HIGH  FREQUENCY  (VHF)  AND  ULTRAHIGH 
FREQUENCY  (UHF)  BAND  SKY  NOISE  ESTIMATES  FROM 
ANALYSIS  OF  34.5  MHz  AND  408  MHz  SKY  SURVEYS 


A.  SKY  SURVEYS 

Several  all-sky  radio  surveys  exist.  Of  these,  two  in  the  frequency  band  of  interest 
have  been  singled  out  for  analysis.  The  first  survey,  at  408  MHz,  was  published  in  1982 
by  Haslam,  et  al.  (Ref.  1).  This  survey,  which  covers  the  entire  celestial  sphere  at  an 
angular  resolution  of  0.85°  x  0.85°,  is  a  compilation  of  data  taken  with  four  radio  tele¬ 
scopes  at  various  locations  on  Earth.  The  second  survey,  at  34.5  MHz,  was  published  in 
1990  by  Dwarakanath  and  Shankar  (Ref.  2).  It  covers  declinations  from  -50°  to  70°. 
Observations  were  made  from  a  single  radio  telescope  in  Gauribidanur,  India. 

Both  data  sets  are  in  the  public  domain  and  can  be  obtained  in  Flexible  Image 
Transport  System  (FITS)  format  from  the  Skyview  Virtual  Observatory  web  site, 
http://skyview.gsfc.nasa.gov.  This  web  site  is  an  astronomical  data  repository  maintained 
by  the  Goddard  Space  Flight  Center  (GSFC). 

Radio  FITS  data  are  expressed  in  Kelvins  (K).  Let  the  measured  brightness  of  a 
radio  source  with  right  ascension  a,  declination  8,  and  frequency  v  be  B.  B  has  units  of 
watts  m"2  Hz"1  rad"2.  The  temperature,  T,  associated  with  this  measurement  is  found 
through  the  low  frequency  limit  of  the  Planck  blackbody  distribution  curve,  namely,  the 
Rayleigh- Jeans  law  (Ref.  3), 

d_2  v2kT  2  kT 

"  2  o2  ’  (U 

C  A 

where  k  is  the  Boltzmann  constant. 

Figures  1  and  2  show  the  temperature  maps  of  the  408  and  34.5  MHz  surveys. 
Temperature  scales  are  given  to  the  right  of  each  image.  Each  image  is  expressed  as  a  rec¬ 
tangular  projection  in  Earth  equatorial  coordinates,  with  the  usual  astronomical  convention 
that  right  ascension  a  increases  to  the  left.  The  size  of  each  image  is  360  x  180  pixels 
(1°  x  1°  resolution).  The  galactic  plane  is  the  sinusoidal  yellow  region  in  each  survey. 
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Figure  1.  All-Sky  Temperature  Survey  at  408  MHz 

Prominent  radio  sources  are  identified  in  the  408  MHz  survey.  The  brightest  dis¬ 
crete  source  in  each  survey  is  the  supernova  remnant  Cassiopeia  A,  which  is  located  at 
(350.29°,  58.55°).  However,  the  brightness  and  extent  of  the  galactic  center  (Sagittarius 
A)  at  (265.52°,  -28.92°)  makes  it  the  dominant  radio  noise  source  in  the  sky. 

The  raw  34.5  MHz  data  must  be  pre-processed,  since  it  contains  values  of  zero 
temperature  outside  the  survey  region  (declinations  above  70°  and  below  -50°)  and  a  few 
isolated  zero  and  negative  temperature  values  within  the  survey  region.  These  values  are 
presumably  dropouts  because  of  data  acquisition  difficulties.  For  convenience,  the  zero 
and  negative  values  have  been  replaced  with  the  spherically  averaged  temperature  over  the 
positive- valued  survey  region  Q.R, 

(T)  =  ^jT(a,8)dQ  ,  (2) 

C1R 
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Galactic  Temperature  (K)  at  34.5  MHz:  Earth  Equatorial  Coordinates 

1 3000000 

1000000 

SSI  3 00000 
1 100000 


30000 


10000 


1000 


60  30  0  330  300 

Right  Ascension 


210  180 


Figure  2.  All-Sky  Temperature  Survey  at  34.5  MHz 

where  dQ-cosSddda.  The  image  shown  in  Figure  2  is  constructed  from  the  pre- 
processed  34.5  MHz  data  set. 


B.  SPECTRAL  INDEX  AND  CONE-AVERAGING 

Note  the  great  difference  in  temperature  scales  between  the  two  surveys.  The  aver¬ 
age  temperatures  for  each  survey  can  be  computed  from  Eq.  (2)  by  allowing  Q,R  to  be  the 
entire  celestial  sphere.  The  averages  are  <Tm>  =  22.13  K  and  <7’34  5>  =  12,317  K.  The 
radio  sky,  therefore,  does  not  behave  as  a  blackbody  whose  temperature  is  constant  for  all 
frequencies.  This  is  because  strong  radio  emitters  are  generally  nonthermal  in  nature. 
They  are  believed  to  be  dominated  by  sources  of  synchrotron  radiation  (Ref.  3). 

A  temperature-frequency  dependence  of  the  form 

T  =  av'r)  ,  (3) 
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where  the  exponent  77  is  known  as  the  temperature  spectral  index  and  a  is  a  proportionality 
constant,  is  often  assumed  to  hold  for  nonthermal  radio  sources.  This  relation  has  a  solid 
theoretical  justification  for  single  synchrotron  radiation  sources  (Ref.  3).  Spectral  indices 
are  typically  in  the  2.4  to  2.7  range. 


For  the  radars  of  interest,  having  an  estimate  of  the  spherically  averaged  sky  tem¬ 
perature  at  frequencies  in  the  30  to  1,000  MHz  range  is  desirable.  This  is  done  as  follows. 
Assuming  that  the  model  [Eq.  (3)]  holds,  the  nominal  spectral  index  77  can  be  found  from 
the  spherically  averaged  temperatures  computed  previously: 


77  = 


ln(<  ?34,5  >  /  <  7^08  >) 

ln(408/34.5) 


=  2.56 


(4) 


The  constant  a  is  then  found  from  Eq.  (3)  by  using  the  average  temperature  at  one 
of  the  two  frequencies  with  77  =  2.56.  The  result  is  a  =  1.06  x  10®  in  units  where  T  is  in 
Kelvins  and  v  is  in  MHz. 


Note  that  other  methods  can  be  used  to  extract  77  and  a.  For  example,  77  can  be 
found  by  computing  a  spectral  index  for  each  pair  of  corresponding  temperature  values 
between  the  two  sky  surveys.  The  mean  value  provides  an  estimate  of  77.  Carrying  this 
out,  one  finds  <77  >  =  2.57,  which  differs  only  slightly  from  the  result  in  Eq.  (4).  The 
parameter  a  can  then  be  found  self-consistently  by  computing  the  mean  of  the  quantity 
a -Tv  <r?>for  the  data  in  both  surveys.  The  result  is  <a>  =  1.76  x  108.  This  result  is  con¬ 
siderably  larger  than  the  value  of  1.06  x  108,  which  was  found  with  the  first  method.  The 
reason  for  this  difference  is  that  the  second  method  treats  all  temperatures  with  equal 
weighting,  while  the  first  method  properly  treats  the  spherical  support  of  the  data.  Since 
our  interest  is  in  average  temperatures  over  the  desired  frequency  range,  we  rely  on  the  first 
method  to  compute  77  and  a. 

It  is  possible  to  estimate  the  standard  deviation  in  the  temperature  at  each  frequency 
and  the  uncertainties  in  the  parameters  77  and  a  by  using  the  following  heuristic  “bootstrap” 
procedure.  First,  the  standard  deviations  at  the  measured  frequencies  are  computed  from 
the  mean  square  temperature 

(r2)  =  -^Jr2(a,5)ufQ  (5) 

and  the  mean  temperature  [Eq.  (2)]  through  the  general  relation 

<j  =  V<r2>-<r>2  .  (6) 

The  results  are  <x34  3  =  19,055  K  and  <r408  =  35.65  K. 
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Suppose  surveys  at  several  frequencies  v,  were  available,  with  corresponding  mean 
temperatures  Tt  and  standard  deviations  or.  Then,  the  parameters  17  and  a  would  be  found 
by  minimizing  the  chi-square  statistic 


a. 


(7) 


This  would  yield  the  best-fit  parameters,  17  and  a,  and  the  uncertainties  in  these 
parameters,  and  <ra.  Knowledge  of  these  uncertainties  would  allow  one  to  compute 
temperature  variances  at  other  frequencies  through  the  propagation  of  variance  equation 


(8) 


assuming  no  cross-correlations  between  rj  and  a.  Computing  the  partial  derivatives  from 
Eq.  (3)  and  dividing  both  sides  by  T 2  leads  to 


^^-KlnvfcrJ 


(9) 


a 


Assuming  this  expression  can  be  used  for  the  two  data  points  under  consideration, 
the  parameter  uncertainties  <Jn  and  oa  can  be  computed  from  the  known  values  of  T  and  cT 
at  34.5  MHz  and  480  MHz.  The  results  are  Ga  =  1.60  x  108  and  C7„  =  0.09. 

A  real  radar  “sees”  a  much  smoother  temperature  profile  than  those  shown  in 
Figures  1  and  2.  Consider  a  receiving  antenna  pointed  directly  away  from  Earth  in  some 
direction  (ac,8c).  The  radar  will  generally  have  a  receiving  pattern  that  subtends  a  solid 
angle,  A Q,  that  is  large  compared  to  the  1°  x  1°  resolution  displayed  in  Figures  1  and  2. 
The  received  energy  is  effectively  an  average  over  the  patch  of  sky  contained  within  AQ 
centered  at  (occ,5c). 

The  temperature  profile  seen  by  a  “typical”  radar  is  simulated  as  follows.  Assume 
the  energy  incident  at  the  radar  antenna  emanates  from  a  patch  of  sky  contained  within  a 
cone  in  the  direction  (ac,5c),  with  cone  angle  /3C  (see  Figure  3).  The  cone-averaged  sky 
temperature  in  the  direction  (ac,5c)  is  given  by  the  average  temperature  over  the  solid  angle 
subtended  by  the  cone: 

T„„(a„6c-,Pc)  =  isjna,5)m-ac,S-Sc-,l3,)d£l  ,  (10) 


where  the  integration  is  carried  out  over  the  entire  celestial  sphere  and  the  kernel  K  is  unity 
for  angles  inside  the  cone  and  is  zero  elsewhere.  Figures  4  and  5  show  the  smoothed  tem¬ 
perature  profiles  at  408  MHz  and  34.5  MHz  for  a  cone  angle  pc  =30°. 
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Effective  Temperature  (K)  at  34.5  MHz:  Earth  Equatorial  Coordinates 


40000 

30000 

25000 

20000 

15000 

12000 

10000 


Right  Ascension 


Figure  5.  Cone-Averaged  Temperature  at  34.5  MHz;  /?„  =  30° 

The  averaging  operation  heightens  the  similarities  between  the  two  spatial  tempera¬ 
ture  distributions.  The  galactic  center  spreads  its  effects  over  a  significant  fraction  of  the 
southern  hemisphere,  while  Cassiopeia  A  stands  out  prominently  at  northern  latitudes. 

As  was  done  previously  for  the  raw  temperature  data,  the  standard  deviation  for  the 
cone-averaged  data  can  be  computed  from  Eqs.  (2),  (5),  and  (6).  In  principle,  the  average 
temperatures  before  and  after  cone-averaging  should  be  the  same  since  the  integral 
[Eq.  (10)]  is  “norm-preserving”;  however,  there  are  slight  changes,  presumably  caused  by 
errors  in  approximating  the  integral  in  Eq.  (10)  with  a  finite  sum.  Table  1  summarizes  the 
averages  and  standard  deviations  before  and  after  cone-averaging.  Cone-averaging  is  seen 
to  reduce  the  standard  deviation  significantly  in  each  survey. 

The  bootstrap  method  discussed  previously  can  also  be  used  to  estimate  the 
parameter  uncertainties  and  the  temperature  standard  deviations  for  the  cone-averaged  data. 
Using  the  data  in  Table  1,  one  finds  on  =  0.09  and  oa  =  9.2  x  107  from  Eq.  (9). 
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Table  1.  Mean  Sky  Temperatures  and  Standard  Deviations 
Before  and  After  Cone-Averaging 


WMM 

'  t  Or 

MSB 

408  MHz  (raw) 

22.13 

35.65 

408  MHz  (cone-averaged) 

22.09 

22.51 

34.5  MHz  (raw) 

12317 

19055 

34.5  MHz  (cone-averaged) 

12295 

11332 

Our  estimate  of  the  sky  temperature  in  the  30  to  1,000  MHz  range  is  taken  to  be  the 
spherically  averaged  mean  temperature  computed  from  Eq.  (3)  plus  one  sigma  computed 
from  Eq.  (9),  with  the  cone-averaged  parameter  uncertainties  <Jn  =  0.09  and  aa  =  9.2  x  107 
(Ref.  4).  Figure  6  shows  this  temperature  along  with  the  mean  temperature  and  the  one 
sigma  level  for  the  raw  data. 


Figure  6.  Estimated  Sky  Temperatures  Over  the  VHF  and  UHF  Bands 
C .  CONCLUSION 

Analysis  of  radio  sky  surveys  at  408  MHz  and  34.5  MHz  has  furnished  estimates 
of  sky  temperature  throughout  the  VHF  and  UHF  bands  (30  to  1,000  MHz).  This  has 
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been  accomplished  through  an  empirical  method  that  extends  the  fitted  powerlaw  tempera¬ 
ture  model  to  the  entire  frequency  range  of  interest.  Estimates  of  noise  variance  have  also 
been  made  using  a  simple  propagation  of  variance  assumption.  Results  of  this  analysis  are 
prerequisites  for  the  evaluation  of  system  noise  temperatures  in  proposed  radars  (Ref.  4). 
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GLOSSARY 


CCC&D 

Counter  Camouflage  Concealment  and  Deception 

DARPA 

Defense  Advanced  Research  Projects  Agency 

FITS 

Flexible  Image  Transport  System 

GSFC 

Goddard  Space  Flight  Center 

IDA 

Institute  for  Defense  Analysis 

K 

Kelvin 

MHz 

megahertz 

UHF 

ultrahigh  frequency 

VHF 

very  high  frequency 

WWW 

World  Wide  Web 
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